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Abstract 

We predict that the combination of space and time modulational instabilities occur- 
ring via parametric wave-mixing in quadratic media leads to coloured conical emission. 
This phenomenon should be observed under conditions usually employed in second- 
harmonic generation experiments. 

Modulational instability (MI) stands for the exponential growth of periodic, long-wave- 
length perturbations occurring either in time (i.e., temporal break-up [1]), or space (fila- 
mentation, e.g. Ref. [2]), due to interplay of Kerr nonlinearity and group- velocity dispersion 
(GVD) or diffraction, respectively. Mis occur also in quadratic media, in the form of para- 
metric instabilities of two-color pump beams, encompassing both the cases of propagation 
invariant solutions (i.e., eigenmodes [3, 4]), and periodic pump evolutions [5, 6] as in, e.g. 
unseeded second-harmonic generation (SHG). Filamentation has been observed both in un- 
seeded and seeded SHG geometries, in 1+2 (bulk) and 1 + 1 (planar waveguides) dimensions 
[4, 5, 7], while temporal MI have not been observed yet. On this ground, one might naturally 
wonder about the coupling of spatial and temporal degrees of freedom, as occurs, e.g., in 
self-focusing of short pulses. In this letter, our purpose is to show that, in the presence of 
both dispersion and diffraction, spatio-temporal MI in quadratic bulk samples cannot be 
simply seen as the mere superposition of filamentation and temporal break-up. Rather, it 
leads to a new phenomenon, namely coloured conical emission (CCE), where rings of differ- 
ent colours are expected to grow exponentially on top of both pump far-fields. As we show 
CCE via SHG has specific features which make it different from that earlier predicted for 
Kerr media [8, 9, 10]. It differs also from other space-time phenomena characteristic of SHG, 
such as snake instabilities [11], where the spatial breakup of SHG solitons associated with 
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temporal MI does not lead to substantial modifications of the beam angular spectra, which 
remain confined nearly on-axis. 

Although the phenomenon of CCE implies component which travel off-axis, the angles 
usually involved are rather small and it is reasonable to use the following paraxial dimen- 
sionless model for SHG in noncritical configuration (i.e., no spatial walk-off) 



9u l . ^1^2 f3 1 d 2 u 1 , 



(i) 



du 2 (7 2 2 du 2 f3 2 d 2 u 2 u\ 

— + -V,u 2 + i 5— - - — + 5ku 2 + — 



Here, the (cw plane-wave) fields are conveniently normalized in terms of their total in- 
tensity If to yield \u 2 \ 2 + |«i| 2 /2 = 1. Accordingly the real-world distance Z = zZ n j is 
fixed by the nonlinear length scale Z n i = (xh)^ 1 , where x = ko(2r]Q/nln 2 ) 1 ^ 2 d e ff [3]. The 
normalized retarded time t = (T — Z/V 9 i)/Tq and the Laplacian = d 2 + d 2 , where 
(x,y) = (X, Y)/Rq, are given in terms of temporal and transverse scales, T = {\k'[\Z n i) x l 2 
and Rq = (Zni/ki) 1 / 2 , respectively. Moreover a m = ki/k m , m = 1,2 (we set a 2 = |), 
Pm = ^m/l^il) = d 2 k/du 2 \ muJ are GVDs, 5k = (k 2 — 2k\)Z n i is the phase mismatch, and 
5 = Z n i/Z w stands for group- velocity mismatch (GVM), Z w = T (l/V g2 — l/V^i) -1 being the 
walk-off length (V gm = dk/dw]^). 

Let us consider propagation-invariant solutions of Eqs. (1), i.e. cw plane- wave eigenmodes 
u m = u^ exp(imi3z), m = 1,2, where uf — f3 — 5k ± (5k 2 + 12) 1 / 2 /6, uf 2 = 2[1 — |iz.^| 2 ] , 
the upper and lower sign standing for in-phase (0 = arg(u 2 — 2u±) = 0, existing for — oo < 
5k < 2) and out-of-phase (0 = ir, oo > 5k > — 2 ) modes, respectively [3]. We study the 
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stability of these waves by inserting in Eqs. (1) the fields 

u m (r, t, z) = [u mQ + a m (r, t, z)] exp(im{3z),m = 1, 2, 
where u m0 = u^o, and we consider a radially symmetric perturbation: 

a m (r,t,z) = /+(r,t)exp(Az) + f~(r, t)* exp(X*z). 



(2) 



(3) 



After linearization, we obtain a system of four partial differential equations in the un- 
knowns /f^O"^)- Without loss of generality, we introduce the Fourier-Bessel transform of 
F^(K, il) = T (r, t)], which permits to express the perturbation amplitudes through 
the inverse transform = T~ x [F^], as a superposition of conical Jo Bessel waves [12] of 
transverse wavevector K = [K 2 + K 2 } 1 ^ 2 and frequency f2, as 



oo poo 



oo JO 



F± (K, Q) J (Kr) e int KdKdtt. 



(4) 



In terms of the vector F = [Ff F 1 F 2 + F 2 ] T of transformed amplitudes, the linearized 
system becomes algebraic and takes the form 



V 



Hi + iX 


U 2 


WlO 







Qi — iX 





UlO 







tt 2 - 5il - 


-iX 










Q 2 + - iX 



F = 0, 



J 



where Q m = /3 m Q 2 /2 — a m K 2 /2 — m(3 + (m — l)5k. The compatibility condition for this 
4x4 system gives the dispersion relation A = X(K, Q). The eigenvalues A have cumbersome 
expressions, except for vanishing GVM (5 = [13]), for which we find (set 2/ = u 20 — 2u\ 

ill ~ H2, fl = (HlH2 ~ uio) 2 ~ l4 ilf), 



,2 

no 
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(5) 



Instabilities sets in when at least one of the four eigenvalues A in Eq. (5) have a positive 
real part, thus causing the perturbation to grow exponentially with gain g = g(Q, K) = 
2Re(A). Two qualitatively different regimes are found depending on the sign of GVD. In 
Fig. 1 we display the typical spectral gain curve obtained in the anomalous GVD regime 
(A, 2 = — 1)- As shown, the gain is relatively narrowband, and its axisymmetric feature (g 
peaks at Q 2 + K 2 = constant) reflects the fact that space and time play indeed the same 
role. Conversely, when the GVD is normal, perturbations with large wavevector (angle) K 
and detuning Q, can be amplified. Three different branches of growing perturbations exist 
even when the collinear process is phase-matched (8k = 0), as shown in Fig. 2(a). This type 
of spatio-temporal MI entails CCE, which can be regarded as the amplification process of 
a superposition of Bessel Jo beams characterized by progressively larger values of of angle- 
frequency (K — Q) pairs. The peak gain [obtained in a wide bandwidth of high frequencies 
K — Q, see Fig.2(a)], is displayed in Fig.2(b) and decreases monotonically for large \8k\ for 
both eigenmodes. 

The underlying physics of CCE involve multiple three-photon processes. Pump photons 
at 2a; travelling on-axis decay into photon pairs at frequency uj — 8uj and uj + 8uj travelling 
off-axis with opposite angles. These in turn can produce off-axis photons at 2uj±8uj via sum- 
frequency mixing with on-axis pump photons at a;. In the limit 8k = ±2 where the pump 
degenerate into the single component uf Q = ±1 (uf = 0), only the downconversion process 
takes place, and we can write the longitudinal projection (the transverse one is identically 
satisfied) of the noncollinear phase-matching condition as k z (2uj)—k z (uj+8uj) — k z (uj—8uj) = 0. 
In the paraxial approximation k z (u ± 80S) = \k(to ± 8uj) 2 — k 2 ] 1 ^ 2 ~ k(u ± Su) — k 2 /2k\ (k t is 
the real- world tranverse wavenumber) and by expanding k(oj ± 8to) = k\ ± Vg{ l 8uj + k'{8u 2 /2, 
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the phase-matching condition yields 

k"5u 2 - k^k 2 = Ah. (6) 

In terms of dimensionless frequencies K = ktRo and Q = 5uTo, Eq. (6) reads as Q 2 — K 2 = 
±2, which turns out to be in perfect agreement with the locus (K — Q) of peak gain as 
obtained from Eq. (5). 

Remarkably we find that CCE is robust against the presence of GVM. In fact for 8^0 
CCE occurs also in the ideal GVD-free limit (flip — 0). The relative magnitude of GVM 
and GVD depends, however, on the intensity I t . For instance, in a LBO crystal pumped at 
A = 1 A*m (k" = 0.045 ps 2 /m, fc£ = 0.25 ps 2 /m) with I t = 50 GW/cm 2 (Z nl = 0.6 mm), we 
find Z w = 0.12 mm which yields 5 = 5.4. Figure 3 shows the different gain branches for this 
case, when the crystal operates in the large mismatch limit k2 — 2ki = —30 cur 1 (8k = —1.8). 
The dimensional gain G = g/Z n i is reported against real-world angles 9 ~ sin# = K/{k\Ro) 
and wavelength detunings AA = AQf2/(27rcTo). 

Though we have considered perturbations growing on top of the SHG eigenmodes, CCE 
can occur also from different launching conditions entailing exchange of energy between u 
and 2u pump beams. This case is similar to spatial or temporal MI occurring in unseeded 
SHG, and can be analyzed by means of Floquet techniques [5, 6]. In this letter, we rather 
present numerical evidence for such phenomenon, by integrating Eqs. (1) in l+l+lD, setting 
d y = 0. As an example Fig.4 shows the spectral content of the two fields after propagation 
(z = 15), as obtained when only the pump beam u\ is launched at z — 0, on top of white 
noise. As shown, CCE favours the amplification of spatio-temporal frequencies which closely 
follow the spectral gain [see Fig. 2(a)]. Due to the large bandwidth of the process we have 
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also found that CCE takes place from short pulse, narrow beam excitations, a case which 
will be discussed elsewhere. 

We end up by emphasizing that CCE in SHG exhibits important differences with respect 
to CCE in true Kerr media. First, SHG makes CCE accessible in materials which are virtually 
GVD-free, and, in the large mismatch limit, for both signs of effective Kerr nonlinearity in 
the same physical material. However, the major importance of CCE stems from the fact that 
it can trigger the formation of novel spatio-temporal propagation-invariant wavepackets of 
Eqs. (1), i.e. u m (r,t,z) = q m (r,t) exp(imj3z) [14]. In the normal GVD regime, q m (r,t) exist 
for (3 < in the form of a polichromatic superposition of Jo beams [12], so-called nonlinear 
X-waves [14], which exist also in the linear limit [15]. In SHG when the Jo perturbations 
are amplified at the expense of an out-of-phase pump beam, they develop the right sign 
of nonlinear phase-shift {(3 < 0) for the formation of a nonlinear X-wave. Viceversa, in 
Kerr media, where nonlinear X-waves also exist, cw plane- waves and localized wavepackets 
experience opposite nonlinear phase-shifts, and such mechanism is unlikely to take place. 

In summary we have shown that spatio-temporal MI of propagation-invariant as well as 
dynamically evolving pump beams in SHG gives rise to conical emission, which might play 
an important role in the space-time dynamics. Such phenomenon is intrinsic to quadratic 
materials and should not be confused with other recently reported mechanisms of conical 
emission [16]. 
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Fig. 1. MI gain g(Q,K) in the anomalous GVD regime for the in-phase eigenmode at 
phase- matching 5 k = 0. 

Fig. 2. Normal GVD regime (Pi,2 — 1) : ( a ) Level plot of conical emission gain g(Q,K) 
for = eigenmode (similar plot holds for <p = 7r); (b) Peak value of g vs. 5k. 

Fig. 3. Level plot of CCE dimensional gain G for the in-phase eigenmode in a LBO 
crystal (see text for data). 

Fig. 4. Spectrally selective amplification of white noise seen at z = 15 from unseeded 
collinearly phase- matched SHG [5k = in Eqs. (1)] pumped by a cw plane- wave pump at 
u: (a) MA',fi)|; (b) \u 2 {K,Q)\. 
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